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Galaxies of different types are not equally distributed in the Local Universe. 
In particular, the supergalactic plane is prominent among the brightest 


ellipticals, but inconspicuous among the brightest disk galaxies. This striking 
difference provides a unique test for our understanding of galaxy and 
structure formation. Here we use the SIBELIUS DARK constrained simulation 
to confront the predictions of the standard Lambda Cold Dark Matter (ACDM) 
model and standard galaxy formation theory with these observations. We 
find that SIBELIUS DARK reproduces the spatial distributions of disks and 
ellipticals and, in particular, the observed excess of massive ellipticals near 
the supergalactic equator. We show that this follows directly from the local 
large-scale structure and from the standard galaxy formation paradigm, 
wherein disk galaxies evolve mostly in isolation, while giant ellipticals 
congregate in the massive clusters that define the supergalactic plane. Rather 
than being anomalous as earlier works have suggested, the distributions 

of giant ellipticals and disks in the Local Universe and in relation tothe 
supergalactic plane are key predictions of the ACDM model. 


The Local Supercluster is the largest structure in the Local Universe’; 
it defines the supergalactic plane and the supergalactic coordinate 
system”. While the plane was originally defined by the relatively nearby 
clusters such as Virgo and Fornax, itis now understood to extend to at 
least redshift z= 0.02 as an excess of bright elliptical and radio galax- 
ies near the supergalactic equator* °. Strikingly, however, there is no 
corresponding excess in bright disk galaxies”*. 

In the standard paradigm of hierarchical structure formation’, 
denser regions are characterized by a greater abundance of dark mat- 
ter haloes and, in particular, more massive ones". Standard galaxy 
formation theory predicts that disk galaxies evolve largely in isolation, 
growing primarily through in-situ star formation fuelled by the continu- 
ous accretion of gas™”. By contrast, the higher ambient temperature 
in denser regions restricts the gas supply into haloes, which leads to 
gas depletion and the eventual quenching of star formation. Denser 


regions are also characterized by a higher merger rate and a higher 
fraction of gas-poor mergers, both of which are understood to lead to 
the formation of elliptical galaxies. 

However, the metamorphosis of galaxies is believed to take 
hundreds of millions of years, precluding its direct observation in 
individual objects. The differing distributions of galaxies of different 
types and masses in different environments within the Local Uni- 
verse thus provide an important opportunity to simultaneously test 
the ability of the cosmological framework to explain the observed 
large-scale structure and to enhance our understanding of galaxy 
formation. 


ALocal Universe constrained simulation 
Direct comparisons of model predictions to observations in the Local 
Universe require a constrained simulation” ’. In this work, we use the 
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Fig. 1| Distributions of disks, intermediates and ellipticals in the Local 
Universe uptoz = 0.02. Hammer projection in supergalactic coordinates of the 
distribution of the most massive disks (left), intermediates (centre) and 
ellipticals (right) in the redshift range of 0.01 < z < 0.02, as reported in the 2MRS 
survey (top) and as modelled in SIBELIUS (bottom). The shaded areas correspond 
to |sin(SGB)| < 0.2 (—11.5° < SGB < 11.5°) around the supergalactic equator. 
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Dark and light symbols show individual galaxies that lie inside and outside this 
region, respectively, and percentages in the top right of each panel express their 
relative numbers. Grey points denote lower mass galaxies (fainter for |b| < 10°) 
irrespective of morphology; labels indicate the positions of five galaxy clusters. 
In both the 2MRS data and the SIBELIUS simulation, a significantly higher fraction 
of massive ellipticals than of massive disks lie close to the supergalactic plane. 


SIBELIUS DARK (SIBELIUS) simulation” that is designed to reproduce 
the galaxy distribution out to a comoving distance of ~200 Mpc (with 
z= 0.04). The initial conditions are constrained using the Bayesian 
Origin Reconstruction from Galaxies (BORG) algorithm” ~ to rep- 
resent the most likely initial conditions that give rise to the observed 
large-scale structure in the 2M++ galaxy redshift survey”, inthe Lambda 
Cold Dark Matter (ACDM) cosmology, with parameters for dark energy 
fraction Q, = 0.693, all-matter fraction Q,, = 0.307, baryonic matter 
fraction Q, = 0.04825, 8 Mpc matter density fluctuations ds = 0.8288, 
primordial density fluctuation spectral index n, = 0.9611 and the Hubble 
Constant Hy = 67.77 km s” Mpc”. 

To predict the SIBELIUS galaxy population, accounting for the 
environment and merger history, the semi-analytical galaxy formation 
model GALFORM”*” was applied to its halo merger tree. In GALFORM, 
galaxy properties such as the stellar mass and morphology are calcu- 
lated based on the properties of the evolving dark matter structures, 
such as mass, accretion rate, spin and mergers, with analytical pre- 
scriptions for processes including gas cooling, star formation, stellar 
evolution, black hole formation, feedback from supernovae and active 
galactic nuclei (AGN), dynamical friction and dynamical instabilities. 
Mergers (particularly major mergers) and starbursts caused by disk 
instabilities are the primary processes that can lead to the growth of 
the elliptical-defining bulge and the transformation of a disk galaxy 
into an elliptical galaxy”. 

Of particular relevance to this study is that the BORG algorithm 
used inthe construction of the constraints is agnostic to the astrophysi- 
cal processes that may give rise to the divergent formation pathways 
of disks and ellipticals, while the GALFORM model is agnostic to the 
particular structures present in the Local Universe, which, however, 
determine the properties of the merger trees. 

We classify the SIBELIUS galaxies according to their bulge-to-disk 
mass ratio, B/D, with distribution as shown in Extended Data Fig. 1. We 
designate those with B/D < 1/4 as disk galaxies, those with 1/4 < B/D < 8 
as intermediates and those with B/D > 8 as ellipticals”°. As shown in 
Extended Data Fig. 2, our results are not sensitive to these particular 
choices. We focus here mostly on the distributions of disks and ellipti- 
cals, which are the subjects of refs. 8,27, but the intermediates provide 
a further point of comparison. 


Predicted and observed distributions 

For comparison to the observations and, in particular, to the results of 
ref. 27, we select from the Two Micron All Sky Survey (2MASS) galaxy 
redshift catalogue” (2MRS) the brightest objects in the redshift range 
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Fig. 2 | Excess of ellipticals relative to disks, as a function of sample size. Top, 
the fraction of the most massive n disks (fp) and the most massive n ellipticals (fẹ) 
in SIBELIUS and 2MRS at 0.01 < z < 0.02, located at |sin(SGB)| < 0.2. Middle, 

Fi:/fp, the fraction of the most massive n ellipticals divided by the fraction of the 
most massive n disks at |sin(SGB)| < 0.2. Bottom, ratio f;/fp, but now for 
|sin(SGB)| < 0.1. As indicated by the red circles, the values of f;/f, in 2MRS peak at 
the value of n chosen in ref. 27. 


0.01 < z< 0.02 by absolute K-band magnitude. Excluding objects with 
absolute galactic latitude |b| < 10° and applying a magnitude limit of 
Mx < -25.19, we obtain an identical sample of galaxies composed of 54 
disks, 53 ellipticals and 73 intermediates. 

In Fig. 1, we compare the distributions in supergalactic coordinates 
as reported in 2MRS, at the top, and as calculated by SIBELIUS, at the 
bottom. From left to right, diagrams show the positions of the most 
massiven = 54 disks (in blue), n = 73 intermediates (in purple) andn = 53 
ellipticals (in red). In all diagrams, grey dots show the positions of 
fainter galaxies tracing the local large-scale structure. Also shown are 
the locations of five galaxy clusters in this redshift range. Orange lines 
enclose a region of supergalactic latitude (SGB) defined by 
|sin(SGB)| < 0.2 (—11.5° < SGB < 11.5°froundthesupergalactic plane. 
Galaxies within and outside this region are denoted by dark and light 
symbols, respectively. 

In both 2MRS and SIBELIUS, the fraction of galaxies close to the 
supergalactic equator depends on their type. This fraction increases 
from disks to intermediates to ellipticals; the most massive ellipticals 
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Fig. 3 | Distributions of SGB in the redshift range 0.01 < z < 0.02 for 
galaxies of different morphology. Disks are shown on the left, intermediates 
in the middle and ellipticals on the right. In each case, galaxies with |b| < 10° 
are excluded. Histograms show the 2MRS data and results from SIBELIUS. The 
top row shows then most massive objects in each class to match the samples 
of ref. 27, the bottom row shows the full 2MRS samples in this redshift range. 
In each case, we compare to an equal number of objects from SIBELIUS. Red 


lines show the isotropic expectation accounting for sky coverage. In both the 
observations and the simulation, the concentration toward the supergalactic 
equator is much more prominent in the most massive ellipticals than in the 
most massive disks; intermediates fall in between. By contrast, for the larger 
sample sizes, the distributions of different galaxy types do not show this effect 
to such an extent. 


are more strongly clustered than the most massive disks. We explore 
the origin of these differences in the next sections. 

Due to the relatively small sample sizes, it is worth considering 
the possibility of a coincidence, or the ‘look elsewhere’ effect. In 
Fig. 2, we compare the fractions of the most massive n disks f, and 
the most massive n ellipticals fz, at |sin(SGB)| < 0.2 in 2MRS and 
SIBELIUS, at 0.01 < z < 0.02, as functions of n. Overall, we find good 
agreement between the simulation and the observations and a 
significant excess f; > fp for all n < 150 (Note that at large n, complete- 
ness may affect the results). We also find that in 2MRS, the relative 
excess peaks at the values used by ref. 27. Comparisons at any other 
n> 25 yield a smaller relative excess, indicating the presence of a 
‘look elsewhere’ effect. In both SIBELIUS and 2MRS, the relative 
excess increases for a narrower band of |sin(SGB)| < 0.1, albeit with 
even smaller sample sizes. 

Drawing random subsamples of 50 out of each of the 100 most 
massive disks and ellipticals in SIBELIUS yields configurations where 
ellipticals are more clustered than disks (that is more objects within 
|sin(SGB)| < 0.2) in approximately 95.2% of cases. This number 
decreases slightly to 93.5% for |sin(SGB)| < 0.25 and to 84.5% for 
|sin(SGB)| < 0.15. For 2MRS, the corresponding numbers are 88.0%, 
86.7% and 87.0%. The excess is equally significant in SIBELIUS and 2MRS. 

In the top row of Fig. 3, we show the distributions of sin(SGB) 
at 0.01 < z < 0.02 for the same numbers of most massive disks (left), 
intermediates (middle) and ellipticals (right) as used by ref. 27. Black 
bars represent the 2MRS data, filled coloured bars represent SIBELIUS, 
for which points denote individual objects. As was already apparent 
from Figs. 1 and 2, we find a strong excess toward low supergalactic 
latitudes of the most massive ellipticals in SIBELIUS, in good agreement 
with 2MRS. There is a much weaker enhancement for the most massive 
disks, also in good agreement with observations. In both cases, the 
distributions of the most massive intermediates fall in between: they 
are more strongly clustered than the most massive disks, but less so 
than the most massive ellipticals. 

When we expand our samples to match the respective sample 
sizes of the 2MRS catalogue, as shown in the bottom row of Fig. 3, 
we find broadly similar distributions for disks, intermediates and 
ellipticals and excellent agreement between the simulation and the 
observations. To assess the significance of the agreement of the SGB 
distributions in SIBELIUS and 2MRS, we generate random permutations 
of the observational data. This gives a probability for the conspicuous 
agreementin the bottom row of Fig. 3 to arise by chance of <0.001% for 


either disks or intermediates and ~0.02% for ellipticals (which have a 
smaller sample size). 

The results of SIBELIUS thus reproduce those presented in ref. 27 
derived from 2MRS: there is a clear difference in clustering between 
ellipticals and disks in relation to the supergalactic plane, but this is 
only manifest among the most massive objects, not inthe much larger 
samples. 


Causes of the separation 

We identify several causes for the different distributions of the 
most massive disks and ellipticals in the SIBELIUS simulation and, 
by implication, in the Local Universe. We first note that the most 
massive ellipticals in SIBELIUS are much more massive than the 
most massive disks, and the most massive galaxies in SIBELIUS are 
overwhelmingly elliptical. For instance, while the 100 most mas- 
sive disks in the range 0.01 < z < 0.02 have a median stellar mass of 
2.9 x 10"° M, (s.d. 0.8 x 10° M,), the corresponding figure for the 100 
most massive ellipticals is 3.0 x 10" M, (s.d. 1.7 x 10" M,), an order of 
magnitude higher. 

Inthe left chart of Fig. 4, we compare the distributions of sin(SGB) 
for galaxies of different stellar mass, irrespective of galaxy type. The 
most massive 10 or 100 galaxies, with stellar masses above 5 x 10" M, 
and 2.5 x 10° M,, respectively, are significantly more concentrated 
toward low SGB than less massive galaxies. 

This mass difference agrees with observational studies of galaxies 
inthe Local Universe”, which also find that the most massive galaxies 
are overwhelmingly elliptical. Using neutral atomic hydrogen and ultra- 
violet data, ref. 32 recently measured the fraction of ellipticals among 
galaxies with stellar masses above 10""* M, to be ~87%. Hydrodynamic 
simulations” and the GALFORM semi-analytical model, both of which 
are calibrated to reproduce the observed mass functions, qualitatively 
reproduce these findings. This suggests that the mass difference we 
identify in SIBELIUS reflects that in the real Universe and explains 
part of the difference in the distributions of sin(SGB), in line with the 
standard model prediction® and observations” that more massive 
galaxies are more strongly clustered. 

The K-band magnitudes in 2MASS for the most massive disks 
and ellipticals are in fact similar”. We attribute this to the previously 
reported systematic underestimation of the magnitudes of massive 
ellipticals in 2MASS**°. Applying the correction suggested by ref. 39 
to the ellipticals would result in a much larger fraction of ellipticals 
among the brightest galaxies, in line with the other studies mentioned 
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Fig. 4 | Distributions of SGB in the redshift range 0.01 < z < 0.02 for 
simulated galaxies and host haloes, ranked according to different criteria. 
Left, distributions of galaxies ranked by stellar mass; the 10 or 100 most massive 
galaxies in the SIBELIUS simulation are significantly more concentrated to low 
supergalactic latitudes compared to less massive galaxies. The distributions of 
the 1,000 or 10,000 most massive galaxies are similar. Middle, distribution of the 
most massive haloes, ranked by Mjo9. The most massive haloes (which contain 


a high number n; of the 50 most massive ellipticals but a much lesser number 

Np, of the 50 most massive disks) are close to the supergalactic plane. Right, 
distributions of ellipticals, disks and intermediates in the mass range of the most 
massive 100 disks, that is, at the same stellar mass. Bars indicate +1o statistical 
error. Even comparing at equal mass, ellipticals and intermediates are more 
strongly concentrated near the supergalactic plane than disk galaxies. 
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Fig. 5| Distributions of galaxies beyond z = 0.02. Hammer projection in 
supergalactic coordinates of the distributions of the most massive disks (blue) 
and ellipticals (red) in 2MRS (top) and SIBELIUS (bottom) in three redshift 
slices. Galaxies are selected at |b] > 10° and with masses similar to the 100 most 
massive galaxies of each type for z in the range 0.01-0.02. As in Fig. 1, dark and 
faint circles denote galaxies within and outside +11.5° of the supergalactic plane, 
respectively. Fainter galaxies are shown in grey. The percentages in the top right 


of each panel express the fraction of galaxies in each class that are within +11.5° 
of the supergalactic plane. The ellipticals more closely trace the distribution of 
the most prominent structures, while the disks are more equally distributed. The 
excess of galaxies at low supergalactic latitudes extends to z = 0.025. At higher 
redshifts, the distributions of both structures and galaxies are no longer aligned 
with the supergalactic plane. 


above and with SIBELIUS. Of course, results based on ranking galaxies 
within the same morphological class remain unchanged. 

We also find evidence that the different environments on or off 
the supergalactic plane act to separate the most massive galaxies by 
morphology, besides simply sorting them by mass. Inthe centre chart 
of Fig. 4, we show the distributions in sin(SGB) of host haloes in different 
mass ranges. As expected, the most massive 10 or 25 haloes in SIBELIUS 
in this redshift range are strongly biased toward the supergalactic 
plane. They also host a significant number n; of the 50 most massive 
ellipticals, but only asmall number n, of the 50 most massive disks. For 
example, 19 of the 50 most massive ellipticals but only 4 of the 50 most 
massive disks are in the 25 most massive haloes. The 50 most massive 
ellipticals reside in 45 haloes with a median mass of Myo. = 9.0 x 10° Ma, 
while the 50 most massive disks are all in different haloes with a median 
mass of Mj, = 1.9 x 10” M,, over 50 times lower. 

In the right chart of Fig. 4, we show the distribution of galaxies of 
the different types—ellipticals, disks and intermediates—in the mass 
range of the 100 most massive disks. All three galaxy samples have very 
similar masses, but the disk galaxies are still significantly less strongly 
clustered than ellipticals or intermediates of similar mass. This further 
emphasizes the role of the environment: galaxies of the same mass are 
less likely to evolve into massive disks in the dense environment that 
marks out the supergalactic plane. 


Inthe Horizon-AGN hydrodynamical simulation, two scenarios for 
the formation of massive disk galaxies are identified”: 1) an exception- 
ally quiet merger history or 2) alate merger of a spheroid with a massive, 
gas-rich satellite; both are less likely to occur in a dense environment, 
where mergers are more common and cold gas fractions are lower. This 
provides an additional explanation, an environmental explanation, 
for the relative paucity of massive disks near the supergalactic plane. 


The extent of the supergalactic plane 
That SIBELIUS is designed to reproduce structures out to ~200 Mpc 
also lets us investigate the extent of the supergalactic plane beyond 
z= 0.02. Figure 5 shows the sky distributions of disks and ellipticals in 
2MRS (top) and SIBELIUS (bottom) in redshift shells out to z = 0.035. 
The overdensity at low supergalactic latitudes only extends to 
z= 0.025. Beyond this, structures no longer align with the supergalactic 
plane, the distributions become less dominated by individual clusters 
and the fraction of galaxies near the equator fluctuates around the 
isotropic expectation value of ~20%. However, in every redshift range, 
both in 2MRS and SIBELIUS, elliptical galaxies are more strongly clus- 
tered than disks. We also investigated the presence of assembly bias, 
thatis, the difference in clustering that is not explained by halo mass“. 
However, as shown in Extended Data Fig. 3, we found no significant 
effect for the comparatively small galaxy samples we consider here. 
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Fig. 6 | Angular correlation functions. The graph shows the angular correlation 
functions (8) of disks (blue) and ellipticals (red) in the SIBELIUS simulation 
(light solid lines) and the 2MRS data (dark dashed lines) in four redshift ranges. 
The top left chart corresponds to Fig. 1, while the top right chart and all bottom 
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uncertainty. In each redshift range within and beyond the Local Supercluster 


and in both 2MRS and SIBELIUS, elliptical galaxies show significantly stronger 
clustering than disks. 


The enhanced clustering of ellipticals compared to disks also 
manifests in the angular correlation functions shown in Fig. 6. In both 
the simulation and the 2MRS data, and in every redshift range, ellip- 
ticals are significantly more clustered than the disks. The peculiar 
excess of massive ellipticals near the supergalactic plane at redshifts 
0.01 < z < 0.02 is not caused by exceptional galaxy formation physics 
but simply reflects the particular structures in this region. 


Discussion 

Constrained simulations such as SIBELIUS provide a unique window 
onto the local matter and galaxy distributions and allow powerful, 
direct tests of cosmological and astrophysical models. While environ- 
mental effects on morphology have been previously demonstrated***°, 
earlier works based on random environments rely on the assump- 
tion that cosmic variance will cancel out. By comparing theory and 
observations in matched environments, we do not have to make this 
assumption. 

The SIBELIUS simulation, in combination with the GALFORM 
semi-analytical model, reproduces the distributions of ellipticals, 
intermediates and disks in the redshift range 0.01 < z < 0.02 and 
beyond. Importantly, it also reproduces the striking difference in the 
clustering of the most massive ellipticals and disks in relation to the 
supergalactic plane. 

We identify two causes for this dichotomy in SIBELIUS. First, the 
most massive ellipticals in the redshift range 0.01 < z< 0.02 are much 
more massive than the most massive disks, and the clustering strength 
depends strongly on galaxy mass. Galaxies which fall between the 
ellipticals and disks in both morphology and stellar mass also follow an 
intermediate spatial distribution. Second, from the fact that the most 
massive disks are less strongly clustered than intermediates or ellipti- 
cals of the same mass, we conclude that the environment prevailing in 
the supergalactic plane inhibits the conditions necessary for massive 
disk formation: a quiet merger history and the continuous supply of 
cold gas. In SIBELIUS, a large fraction of the most massive ellipticals are 


close to the supergalactic plane by virtue of residing within the galaxy 
clusters that define it. 

The GALFORM galaxy formation model is idealized and cannot 
capture the full complexity of astrophysical processes at play in hydro- 
dynamic simulations, let alone in the real universe. Nevertheless, it is 
patently sufficient to reproduce the distinct populations of elliptical 
and disk galaxies observed in the Local Universe. 

The strikingly different distributions of bright ellipticals and disks 
in relation to the supergalactic plane do not require physics beyond 
the standard model. They arise naturally inthe ACDM framework and 
the standard model of galaxy formation represented by the SIBELIUS 
simulation of the Local Universe. They area direct consequence of the 
differences in mass and in the associated bias and of the transformation 
of disks into ellipticals within the clusters that define the supergalactic 
plane. Rather than as an anomaly, the observed distributions emerge 
as a prediction of the ACDM paradigm and as an important benchmark 
for any alternatives. 


Methods 

The SIBELIUS DARK simulation 

Our work is based on the SIBELIUS simulation, a collisionless cos- 
mological simulation which embeds a high-resolution constrained 
region of radius -200 Mpc centred on the Local Group within a1 Gpc? 
periodic volume. SIBELIUS is set up with cosmological parameters 
A, = 0.693, O,, = 0.307, Q, = 0.042825, 0, = 0.8288, n, = 0.9611 and 
Hy = 67.77 km s* Mpc”. The high-resolution region is sampled with 
approximately 130 billion dark matter particles of mass 1.15 x 10’ M, 
with amaximum physical softening length of 1.4 kpc. More details of 
the numerical setup are given in ref. 19, which also includes additional 
comparisons with observations. 

The initial conditions for SIBELIUS are created using white noise 
fields obtained through the BORG algorithm”° ”, which takes as an 
input the observed three-dimensional density field, inferred from the 
2M++ galaxy survey”. BORG is designed to infer the most likely initial 
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configuration that gave rise to the input galaxy density distribution, 
subject to the specified cosmology and an assumed galaxy bias model; 
it represents the most likely underlying matter distribution and forma- 
tion history of the structures in the Local Universe. It is important to 
note, however, that the input galaxy distribution used in constructing 
the initial conditions for SIBELIUS does not differentiate by galaxy mor- 
phology or colour. The reproduction of the observed local large-scale 
structure is thus ‘by construction’, but the reproduction of distinct 
galaxy subpopulations is not. 

Outputs of SIBELIUS are stored 200 times between z= 25 
and z= 0. For each of these snapshots, haloes and subhaloes are 
respectively identified using a Friend-of-Friend algorithm and the 
Hierarchical-Bound-Tracing algorithm”. Subhaloes are further pro- 
cessed into a merger tree, to which the GALFORM semi-analytical model 
is applied, as discussed below. 


The GALFORM semi-analytical model 

SIBELIUS only traces the evolution of matter via gravity. To derive a 
prediction for the galaxy population that can be compared to observa- 
tions, the separate evolution of baryons and the formation of galaxies is 
followed using the semi-analytical model GALFORM. All semi-analytical 
models make the same basic assumptions: that baryons follow the 
accretion of matter (including its angular momentum) onto collapsed 
dark matter haloes** and that the evolution of baryons inside dark mat- 
ter haloes and during mergers and interactions can be described bya 
set of coupled, non-linear differential equations”***°. In GALFORM, 
these equations describe and account for various astrophysical pro- 
cesses including gas cooling, dynamical friction, star formation, stellar 
evolution, black hole formation and evolution and the formation and 
feedback from AGN, stars and supernovae. 

The parameters of the model are set either from ab initio considera- 
tions, or via calibration to observations of the galaxy population. Since its 
first introduction” and through several subsequent updates", when 
applied to ACDM merger trees, GALFORM has been shown to reproduce 
many of the observed properties of the galaxy population, such as the 
neutral atomic hydrogen or stellar mass functions, the ratio of ellipticals 
to disks and the Tully-Fisher relation. It has also been shown to reproduce 
the galaxy population and the properties of individual galaxies (albeit 
with considerable scatter) in the EAGLE hydrodynamic simulation”. The 
version of GALFORM used for SIBELIUS” incorporates different initial 
stellar mass functions during starbursts, feedback that suppresses gas 
cooling in massive haloes and anew empirical star formation law in gal- 
axy disks based on molecular gas content. While it still lacks some of the 
complexity of a full hydrodynamic simulation, GALFORM implements 
our current understanding of galaxy formation physics. 

Importantly, the calibration of its free parameters has been per- 
formed using merger trees from randomized initial conditions with- 
out consideration of galaxy clustering” and independently of the 
constraints used in SIBELIUS, or of the particular structures present 
inthe Local Universe. 


Separation of disks, ellipticals and intermediates. As shown in 
Extended Data Fig. 1, the GALFORM semi-analytical model results in 
a bimodal distribution of galaxy morphologies as characterized by 
the bulge-to-disk mass ratio B/D, or equivalently, bulge-to-total mass 
ratio B/T. We use B/D to classify SIBELIUS galaxies into disks, ellipticals 
and intermediates. By default, we designate galaxies with B/D < 1/4 as 
disks, those with 1/4 < B/D < 8 as intermediates and those with B/D > 8 
as ellipticals”. These two limits are indicated by dotted vertical linesin 
Extended Data Fig. 1. Using these values, 90% of disks have B/D < 0.13, 
while 90% of ellipticals have B/D > 12.6, indicative of a strongly bimodal 
distribution. It can also be seen that the morphology distribution 
depends on mass. In agreement with observations”, the lower mass 
galaxy population contains a larger fraction of disks, while higher mass 
galaxies are more likely to be elliptical. 


Because the thresholds separating the three classes are still some- 
what arbitrary, we also investigate the effect of choosing different 
values. In Extended Data Fig. 2, we repeat Fig. 2 using B/D < 1/A for disks 
and B/D > A for ellipticals, for values of A = 3, 4, 5and 6. The behaviour, 
in terms of an excess clustering of the most massive ellipticals relative 
to disks, is similar in every case, demonstrating that our results are not 
sensitive to the precise delineation between galaxy types. 


The Two Micron All Sky Survey redshift catalogue 

The observed distributions of galaxies are drawn from the 2MRS** 
catalogue, which lists 44,599 galaxies. We include in our analysis all 
galaxies for which redshift and morphological types are available. 
Galaxy morphologies are as assigned in ref. 28, and we follow Peebles’ 
classification, whereby types 21 are classified as ‘spirals’ or ‘disks’, and 
types <-5 are classified as ‘ellipticals’. We additionally identify inter- 
mediate types as ‘intermediates’. 

To sort galaxies by luminosity, we infer absolute K-band magni- 
tudes using the distance modulus. We adopt H, = 70 km s* Mpc” for 
consistency with ref. 27, but we obtain the same galaxy samples witha 
value of Ho = 67.77 kms“ Mpc ‘as used in the SIBELIUS simulation. We 
also adopt the same galactic latitude threshold |b| > 10° to obtain, from 
the 180 brightest galaxies in the redshift range 0.01 < z < 0.02, the same 
samples of 54 brightest disks and 53 brightest ellipticals. We classify the 
remaining 73 as intermediates. We use the same morphology definitions 
and galactic latitude thresholds also at higher redshifts, up toz = 0.035, 
and we apply the same galactic latitude cuts to the simulation data. 

To perform the comparison in redshift space, we use the measured 
redshifts for the 2MRS sample and use a value of H, = 67.77 km s* Mpc? 
(in addition to the proper motions) when assigning redshifts to the 
galaxies in SIBELIUS based on their comoving simulation coordinates. 


Magnitude bias for ellipticals. As noted by Peebles, the absolute 
K-band magnitudes of the brightest 53 disk galaxies and the bright- 
est 54 elliptical galaxies at 0.01 < z < 0.02 in 2MRS are similar, which 
would indicate similar stellar masses”. This is in contrast to other 
observations, to the results of hydrodynamic simulations and also to 
the results of the GALFORM semi-analytical model, according to which 
the most massive ellipticals are significantly more massive than the 
most massive disks and, by implication, the most massive galaxies are 
overwhelmingly elliptical. We attribute this difference to the known 
fact that 2MASS systematically and significantly underestimates the 
luminosity of bright bulges and ellipticals**“°. 

In particular, by comparison to deeper K-band observations of the 
same objects”, the finding was that 2MASS underestimates the luminosity 
of ellipticals and bulges by an average of 0.34 magnitudes, with the effect 
increasing for bright ellipticals. Applying a (somewhat conservative, but 
also crude) correction of -0.34 tothe magnitudes ofall ellipticals changes 
the distribution of morphological types among the brightest galaxies. 
Using the same selection as ref. 27, among the 180 brightest galaxies in 
2MRS with 0.01 < z< 0.02 and |b| > 10°, the number of ellipticals increases 
from 53 to 89, while the number of disks decreases from 54 to 34. Among 
the brightest 100 galaxies, the number of ellipticals increases from 36 to 
56, while the number of disks decreases from 25 to 14. 

It is important to note that this correction would not change the 
ordering of galaxies within the same class; the identity and positions 
of the brightest n ellipticals and the brightest n disks remain the same. 
We have no need to apply a magnitude correction in this work, but its 
effect may explain why, when the 2MASS K-band magnitudes were 
taken at face value, stellar mass and the associated formation bias was 
excluded as a factor for explaining the different distributions of bright 
ellipticals and disks”. 


Angular correlation functions 
As noted in the main text, the strikingly different distributions of bright 
ellipticals and disks in relation to the supergalactic plane reflects the 
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fact that, inthe redshift range 0.01 < z < 0.02, the Local Supercluster is 
the most prominent feature and aligns with the supergalactic equator. 
In other redshift ranges, we find no strong excess of ellipticals near the 
supergalactic plane, but as shown in Fig. 6, the difference between the 
correlation functions of disks and ellipticals persists at all redshifts, 
bothin 2MRS and in SIBELIUS. The astrophysical processes that lead 
to different clustering of disks and ellipticals are universal, but the 
prominent excess of ellipticals at low supergalactic latitudes in the 
redshift range 0.01 < z < 0.02 is caused by the presence of the Local 
Supercluster in this region of the Local Universe. 

To compute the angular correlation functions w(@), we use the 
Landy-Szalay estimator”. For each redshift range and each galaxy 
type, we use the full 2MRS sample, choosing matching numbers of 
most massive galaxies in SIBELIUS and of random points. We apply the 
same |b| > 10° mask for each data set and for the random samples. We 
use the jackknife method to account for the random nature of galaxy 
formation in every data set and also repeat the random sampling to 
estimate the statistical uncertainty. 


Assembly bias 

Elliptical galaxies populate much more massive haloes than disk galax- 
ies in our simulation. As is evident from Figs. 3 and 4, respectively, the 
stronger clustering of ellipticals toward the supergalactic plane for 
0.01<z< 0.02 is present both when classifying galaxies by morphology 
and by host halo mass. 

This may raise the question of whether the stronger clustering of 
ellipticals is driven entirely by the mass of their host haloes, or whether 
there is an additional effect, known as assembly bias“. Assembly bias 
can be measured by comparing the correlation functions of samples 
of galaxies separated by morphology to those obtained when the 
position of each galaxy in the sample is randomly assigned to a halo 
of similar mass, regardless of morphology. A systematic difference 
between the correlation function of the original sample and that of 
the shuffled counterpart must be attributable to factors other than 
halo mass, such as environment or assembly history; such difference is 
referred to as assembly bias. Ref. 41 found that assembly bias accounts 
for approximately 10% of the total bias, being positive for ellipticals, 
slightly negative for disks and decreasing in magnitude for the most 
massive galaxies. 

We performed a similar test for the galaxy populations in SIBEL- 
IUS, showing the resulting angular correlation functions for different 
sample sizes of disks and ellipticals in Extended Data Fig. 3, forthe same 
redshift ranges as in Fig. 6. When comparing to Fig. 6, it is important 
to note that because only central galaxies are considered, the correla- 
tion of galaxies within the same halo (the so-called one-halo term) is 
suppressed, resulting in weaker overall correlations. 

In every case, as expected, we find a stronger correlation 
for the central ellipticals than for the central disks. There is some 
sign of positive assembly bias for our largest samples of galaxies 
(n=2,000), but for smaller samples, the assembly biases can appear 
as both positive and negative and do not appear to be statistically 
significant. 

We have no reason to believe that assembly bias does not also 
apply in the Local Universe, but it appears that the standard halo bias 
is the most significant effect in explaining the different distributions 
of the most massive disks and ellipticals with respect to the superga- 
lactic plane. 


Data availability 

All data used in this work is publicly available. The 2MRS catalogue 
is provided in ref. 28. The halo and galaxy catalogues of the SIBEL- 
IUS simulation can be accessed via SQL at: https://virgodb.dur.ac.uk. 
SQL scripts to obtain the data used in this work are provided with our 
analysis code. The data is also available at https://doi.org/10.6084/ 
m9.figshare.22250602. 


Code availability 

The analysis in this paper was performed in python3 and makes exten- 
sive use of open-source libraries, including Matplotlib v.3.7.0 (ref. 57), 
Scipy v.1.9.3 (ref. 58), NumPy v.1.23.5 (ref. 59), Astropy v.5.1.1 (ref. 60) 
and Corrfunc v.2.5.0 (ref. 61). A documented Juypter notebook con- 
taining the code to reproduce all figures in this paper is available at: 
https://github.com/TillSawala/supergalactic. 
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Extended Data Fig. 1| Distribution of galaxy morphologies as parameterized 
by B/D or B/T in the Sibelius simulation in the redshift range z = 0.01 - 0.04. 
The blue curve shows galaxies in the stellar mass range 10° - 10"°M.,, the red 
curve shows galaxies above 10'°M,,, the grey dashed curve shows all galaxies 
above 10°M,,. There is a clear bimodality in the bulge fraction, and the fraction of 
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galaxies with high bulge fraction increases with increasing stellar mass. The two 
vertical dotted lines denote B/D = 1/4 (B/T=1/5), our default threshold separating 


disks and intermediates, and B/D = 8 (B/T= 8/9), our default threshold separating 
intermediates and ellipticals. 
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bulge-to-disk ratios, less than 1/A for disks and greater than A for ellipticals. 
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type. Note that the correlations are weaker than those in Fig. 6 due to the absence 
of the one-halo term. Assembly bias manifests itself as a difference between the 
original and shuffled samples, but we detect no clear effect of assembly bias for 


these sample sizes. 
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